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Abstract — A system reconfiguration problem is considered 
for three-phase smart power distribution networks featuring 
distributed generation, with the objective of finding the topology 
that minimizes the overall active power loss. In lieu of binary line 
selection variables, the notion of group sparsity is advocated to 
re-formulate the nonconvex distribution system reconfiguration 
(DSR) problem into a convex one. Using the duality theory, it is 
shown that the line selection task boils down to a shrinkage and 
thresholding operation on the line currents. Further, numerical 
tests illustrate the ability of the proposed scheme to identify the 
optimally meshed, weakly-meshed, or even radial configurations 
by adjusting a sparsity-tuning parameter in the DSR cost. 
Constraints on the voltages are investigated, and incorporated 
in the novel DSR problem to effect voltage regulation. 

Index Terms — Smart distribution systems, microgrid, system 
reconfiguration, convex programming, compressive sampling. 

I. Introduction 

The fundamental objective of distribution system recon- 
figuration (DSR) schemes is to identify the topology of a 
distribution network that minimizes the overall active power 
loss 0, 0, 0. DSR byproducts include balancing the 
network load [4|, increasing the system security, and prompt 
(possibly network-wide) power delivery restoration in case of 
localized network failures. Computationally-affordable DSR 
schemes are increasingly advocated in modern distribution net- 
works to enhance their efficiency and stability in the presence 
of distributed generation (DG), energy storage devices, as well 
as dispatchable and elastic loads 0. 

Changes in the network topology are effected by opening or 
closing tie and sectionalizing line switches. These switching 
operations can be either performed manually in situ, or, com- 
manded remotely by a network controller. Thus, the DSR task 
is traditionally approached by associating a binary selection 
variable with each switch @, 0, 0, 0, (TO), (TTJ. Unfor- 
tunately, this choice renders the resultant topology selection 
problem NP-hard lfl2ll . and thus challenging to solve optimally 
and efficiently. This explains why heuristic schemes have been 
largely employed to select the status of the switches. For ex- 
ample, all switches are initially assumed closed in e.g., 0, 0, 
and then they are progressively opened to eliminate (part of) 
the loops. A search over relevant radial configurations based 

This work was supported by the Inst, of Renewable Energy and the 
Environment (IREE) grant no. RL-0010-13, Univ. of Minnesota. Authors are 
with the Digital Technology Center and the Dept. of ECE, University of 
Minnesota, 200 Union Street SE, Minneapolis, MN 55455, USA. E-mails: 

{emiliano, georgios}@umn . edu 

This work has been submitted for publication to the IEEE Transactions on 
Power Delivery on March 20, 2013. 



on approximate power flow methods is advocated in [2 |. Alter- 
native methods rely on fuzzy multi -objective [ 13 1, branch-and- 
bound [6], and genetic algorithms [9|. Off-the-shelf solvers for 
mixed-integer nonlinear programs were employed in e.g. ifTOll 
to solve a joint DSR and optimal power flow (OPF) for 
balanced systems, while Newton methods in conjunction with 
branch-selection heuristic techniques were employed in 0. 
However, computational complexity of these schemes may 
become prohibitive as the size of the network increases, and no 
optimality guarantees are in general available for the obtained 
configuration. 

The present paper leverages contemporary compressive 
sampling tools lfT4l . lfT31 to bypass binary optimization vari- 
ables, and formulate a novel convex DSR problem. Convexity 
here implies that interior point solvers can be utilized to 
solve the challenging DSR problem efficiently and optimally 
even for distribution networks of large size, and with densely 
deployed line switches. DSR solvers able to find a new 
configuration in a few seconds (or even less that one second, as 
shown for the IEEE 37-node feeder |[T6l ) are instrumental for 
network operators to quickly re-shape the distribution grid in 
case of localized outages ifTTIl . and to gauge optimal topologies 
in case of abrupt load or generation variations. 

The proposed convex formulation hinges on the notion 
of group- sparsity, an underlying attribute of the currents 
flowing on the phases of distribution lines equipped with 
switches. This group-sparse problem structure allows one to 
discard binary optimization variables, and select the states 
of the switches by augmenting the DSR cost with a convex 
sparsity -promoting regularization function 1141 . |fT31 . As in 
conventional (group) sparse linear regression, it is shown here 
that the line selection task boils down to a shrinkage and 
thresholding operation on the line currents. This is further 
corroborated through experiments on a modified version IEEE 
37-node feeder lfl6l . where a meshed, weakly-meshed, or 
radial configuration is obtained by simply adjusting a sparsity- 
tuning parameter. Differently from DSR approaches applicable 
to balanced distribution networks 0, 0, fllj], ifTTTl . the 
problem formulated here accounts for unbalanced loads and 
non-zero off-diagonal entries of line admittance matrices. 

Unfortunately, PQ loads and DG units (modeled as PQ loads 
as well) involve challenging nonlinear power flow relations. 
However, since the aim here is to develop a DSR scheme that is 
computationally efficient and yet able to reliably discard lines 
involving high active power losses, the powerful approximate 
load model developed in ifTSl is advocated, and tailored to 
the three-phase setup. Although this load model introduces an 
approximation error (which becomes negligible for sufficiently 
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large values the nominal voltage as shown analytically in 1 18 1, 
and further corroborated numerically here), the payoff here is 
huge, since a convex DSR problem can be formulated even in 
the presence of PQ loads. 

II. Preliminaries and problem formulation 

Consider a portion of the power distribution grid located 
downstream of the distribution substation, that supplies a 
number of industrial and residential loads, and may include 
DGQ The considered three-phase network is modeled also as 
a directed graph (Af,£), where N nodes are collected in the 
set Af := {1, . . . , N}, and overhead or underground lines are 
represented by the set of (directed) edges £ := {(m, n)} C 
Af x Af. Let node 1 represent the point of common coupling 
(PCC), taken to be the distribution substation. Distribution 
systems typically have tie and sectionalizing switches, whose 
states (closed or open) determine the topology of the network. 
Thus, let £ r C £ collect the branches equipped with control- 
lable switches. 

Let V„ m Q {<W, b mn , c mn } and V n C {a n , b n , c n } denote 
the set of phases of line (to, n) 6 £ and node n G Af, 
respectively; I mn G C the complex current flowing from node 
m to node n on phase <j>\ I m G C the current injected at node 
m € Af and phase <fi G V m \ and, V£ G C the complex line-to- 
ground voltage at the same node and phase. Lines (to, n) € £ 
are modeled as n -equivalent components |fl9l Ch. 6], and 
the \P m n\ x \P m n\ phase impedance matrix is denoted by 
Z mn G C' 73 " 1 ™'*' 73 ™"'. Matrix Z mn is symmetric, full-rank, 
and it is obtained from the line primitive impedance matrix 
via Kron reduction lfl9l Ch. 4]. Using Z mn , the \V m n\ x 1 
vector i mn := [{Im n } c/>ev mn ] T collecting the currents flowing 
on each phase of line (m, n) G £ can be expressed as 



Kirchhoff's current law can be written per node n as 



imn — Z m „ {[Vm\Vmn ~ [ v n]p m „) 



(1) 



with v m := [{V„i}^ r J T . Let i„ := [{^} 0eP J r be 
the vector collecting the currents injected at node n, and 
{e^]-0 e -p n and {e^ n } fieVmn me canonical bases of W v ^ and 
W^"*"', respectively. Further, per-node n G Af, define the 
\V n \ x \V mn \ matrix A m l „ := J2 f€Vn I WeTW e*(ef„„) r . 
Suppose for brevity, that the entries of the line shunt ad- 
mittance matrix are negligible (they are, in fact, typically on 
the order of 10 — 100 micro Siemens per mile [16]). On the 
other hand, consider at the expense of minimally increasing 
complexity, that perceptible effects of shunt admittances can 
be readily accounted for in the ensuing problem formulations. 
Under these conditions, and using the definition of A^-™ , 

' Notation: Upper (lower) boldface letters will be used for matrices (column 
vectors); for transposition; (■)* complex-conjugate; and, complex- 
conjugate transposition; denotes the real part, and the imaginary 
part; j = represents the imaginary unit; and Ir.\ is the indicator 
function = 1 if x is true, and zero otherwise). \V\ denotes the 
cardinality of set V; and represent the space of the Af X 1 real 
and complex vectors, respectively. Given a vector v and a matrix V, [v]-p 
denotes a \V\ X 1 sub-vector containing entries of v indexed by the set V, 
and [V]-p, j> 2 the \P\\ X \"Pi\ sub-matrix with row and column indexes 
described by V\ and Vq.. Further, ||v||2 '■= V v^v denotes the £2 norm of 
v. Finally, OjyfxiV and ImxK denote M X N matrices with all zeroes and 
ones, respectively. 



, A (n) i- 



r 



E 



An) 
Nifc 



Ink 







(2) 



where 7V_j. n := {j : (j, n) G £} and Af n ^ := {k : (n,k) G 
£}, respectively. Clearly, i„ = 0\-p n \xi if neither loads nor 
DG units are connected at node n. 



Let S„ 



— \u 



'mnWPmn denote the total complex power 
injected on line (to, n) from node m. If no power is dispelled 
through the distribution line (to, n) G £, S mn coincides with 
the total power transferred to node n; that is, S mn — —S nm . 
A necessary condition for this identity to hold is to have 
an identically zero line impedance matrix. In fact, it readily 
follows from ^ that S mn + S nm = i^„Z mn i m „, and thus 
the total active power loss on line (to, n) G £ amounts to 

+ 9 r {i„ m }3fi{Z m „}3{i m „}. (3) 

However, since typical values of Z mn in overhead and un- 
derground distribution segments render AP mn not negligi- 
ble ITT61 . |fl9ll , it is desirable to select the topology (meaning 
the states of switches on lines £ r), and adjust the complex line 
currents traversing the selected lines, so that the overall 

real power loss Yl( m n)e£ APmn is minimized Q), 0- 

Similar to various DSR renditions J6), Q, H, (9), ifTOl , 
ifTTTl , the topology selection problem will be first formulated 
in the ensuing subsection using binary line selection variables. 
However, since lines may be non-transposed and the spacings 
between conductors may be non-equilateral [19], the off- 
diagonal elements of Z mn are non-zero |[T6l . |[T9l . Thus, 
different from DSR schemes tailored to balanced distribution 
networks (as in e.g., [2|, [8|, IfTOl . ifTTIl '). the problem formu- 
lated here is able to capture current-coupling effects on the 
distribution lines. 

A. Plain-vanilla DSR formulation 

Suppose for the moment that loads are modeled using ideal 
current generators (that absorb current from the network). PQ 
loads and DG power injections (which follow a constant PQ 



(PI) min 

{w},Ma 
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model as well [20 1) will be considered in Section MT-B] Similar 
to (6), 0, H, 13, HO), EU, associate a binary variable 
£77771 G {0,1} with line (m,n) G £r, and suppose that this 
distribution segment is used to deliver power to the loads if 
Xmn = L In this case, the DSR problem can be formulated 
as follows [cf. |2]i] 

^ft {i mn }?R{Z mn }!IR{i mn } 

(m,n)e£ 

+ S r {i m „}3?{Z mn }3{i mn } (4a) 

subject to Q, and 

K 2 {[W]4 + ^ 2 {[iL]4 < C77 , (m, n) G £\£ R (4b) 
^ 2 {[imnU} + 3 2 {[i m „W < C a „ x ^77777 , (m,n) G £ R (4c) 

G {0, 1} , (to, n) G £r (4d) 

where constraints (|4b]i-(|4c|) are enforced on each phase <f> G 
Vmn, and > is a cap for \I mn \ 2 . Clearly, when 
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Xmn — 0, i mn is forced to zero by ( |4c| ), thus implying that 
line (m,n) € £ fl is not used (see also (6), Q, flTUl , ifTTTl ). 
When the desired topology is radial, additional constraints are 
present ifTOll . ITTl . In particular, suppose that the graph (AT, £) 
contains Nj cycles, and collect in the set C L the lines forming 
cycle i = 1, . . . , Nj. Then, to obtain a tree network, it suffices 
to add in (PI) the constraint J2( m ,n)eCi x mn < — 1 per 
cycle i. 

Matrices {$t{Z mn } G M} Vmn I x \ Vmn I } are typically positive 
definite and full-rank (see e.g., the real test cases in lfl6l ): thus, 
it follows that the DSR cost pa} is strictly convex. However, 
presence of the binary variables {x mn } renders (PI) a mixed- 
integer quadratic program (MIQP), which is nonconvex and 
NP-hard [12|. Finding its global optimum requires solving a 
number of subproblems (one per switch status) that increases 
exponentially (2l £n l) in the number of switches. This explains 
why heuristic schemes have been largely employed to select 
the network topology 0, 0, H, US. Alternatively, off-the- 
shelf solvers for mixed-integer (non)linear programs [ 1 1 and 
genetic algorithms [9| have been employed. However, since 
these solvers are in general computationally-heavy, they are 
not suited for optimizing the operation of medium- and large- 
size distribution networks, or, for finding a post-outage system 
configuration in order to efficiently restore power delivery 
network- wide ifTTl . 



In the ensuing Section III compressive sampling tools will 



be advocated to bypass binary selection variables, and re- 
formulate the DSR problem into a convex one. But first, an 
approximate yet powerful load model is outlined next. 



B. Approximate load model 

Consider the well-established "exponential model" relating 
injected (or supplied) powers with voltages {V 7 f } and currents 
{/*} (see e.g., ED) 



V N 



eV n ,neN\{l} (5) 



where Vn is the nominal line-line voltage magnitude of the 
system (e.g., 4.8 kV for the IEEE 37-node feeder [16)); S* is 
the net complex power that would be injected on phase (f> of 
node n if V$ were equal to the nominal voltage -^e^", 

<pff G {0°, -120°, 120°}; and, K n G {0,1,2} is a model 
parameter. Specifically, constant PQ, constant current, and 
constant impedance loads are obtained by setting n n = 0, 
K n = 1, and K n = 2, respectively. 

DG units are typically modeled as constant PQ loads (that 
supply power); thus, let S C N be the (sub)set of nodes 
where DG units are present, and sa,n : = [{Sq n }] T the vector 
collecting the complex power S G supplied by these units on 
each phase of node n. Notice that multiple DG units may be 
present at each node; if this is the case, S G can be readily 
replaced by X)«=i S G n u , with Nu the number of DG units at 
node n. Per phase cf> <G V n , let denote the complex power 
demanded by a wye-connected load at the bus n. Finally, de- 
fine sl,ti := [{S^ n }] T , and suppose as usual that the voltages 



at the substation Vi := [#e j0 ° , #e^ 120 ° , #eJ 120 °l r are 
taken as reference for the phasorial representation 1191 . 

As elaborated further in Section [IV] the use of the nonlinear 
relation |5]l in (PI) introduces an additional source of non- 
convexity [ 10 ], and would render the DSR formulated in the 
next section nonconvex. This, in turn, would exacerbate the 
problem complexity, and would make optimality claims on the 
obtained topology difficult to establish. Instead, the aim here 
is to develop a DSR scheme that is computationally efficient 
yet able to reliably (and optimally) discard lines involving high 
active power losses. Developing such a scheme is instrumental 
for network operators to quickly re-shape the distribution grid 
in case of e.g. abrupt load or generation variations, and to 
promptly restore the power delivery network-wide after an 
outage event ifTTll . 

To this end, the powerful approximate load model derived 
in ifTHIl is advocated to relate injected currents to complex 
powers linearly. Specifically, upon defining the vectors t„ := 
[n T {i n }^ T {i n }} T and (7 Q n := [3? r {s G ,„}, 3 r {s G ,„}^], 
it follows from [18| that the current injected at node n can be 
approximated as 



73 
V N 



-3f?{*„} 



G,n 



(6) 



with <fr„ := diagde^wj-^g-p^), and trg n = for n G 
A/"\({1} U S). The approximation error incurred by |6| is 
infinitesimal for large nominal voltages Vn as analytically 
shown in |18|, and further corroborated by the tests on the 
IEEE 37-node feeder reported in the Appendix. Specifically, 
the average error is just on the order of 0.1 Ampere (a 
relative error of less than 2% - which yields an error on the 
powers that is on the order of the load prediction error). The 
motivation behind |6| is twofold: first, using |6]) in conjunction 
with compressive sampling methods, offers the advantage of 
a convex DSR problem for a system featuring PQ loads. 
Further, from an DSR standpoint, this approximation does not 
jeopardize the ability of the methods proposed in the ensuing 
sections to efficiently capture the topology yielding the lowest 
power loss. Based on this topology, voltages and currents 
are as usual fine-tuned in a subsequent optimization stage by 
employing more sophisticated techniques such as OPF 

III. DSR VIA GROUP-SPARSITY 

Collect first the real and imaginary parts of i mn in the 
vector £ mn := [5R r {i m „}, 3 r {i m „}] r G K 2 I^I, and define 



matrices Z n 



h ® K{Z m „}, A 



(■„.) 



A^n, and 



Mt n := I 2 ® et n (et n ) T , where I 2 is the 2 x 2 identity 
matrix. Key to obtaining a convex re-formulation of (PI) is to 
notice that the entries of £ mn are all zero if line (to, n) G £r 
is not used to deliver power to the loads. In compressive sam- 
pling, this translates to having the vectors {£ mn \(m, n) G £r} 
being group-sparse flU ; meaning that, either £ mn = 0, or, 
the elements of £ mn may all be nonzero. One powerful way 
to capitalize on this attribute of currents flowing on lines 
equipped with switches, consists in augmenting the cost (Hak 



4 



with the following sparsity -promoting regularization term [14] 

r({$mn})-= E X Umnh (7) 

(m,n)gf_R 

where A is a real positive constant. Then, using |7]), and 
discarding binary line selection variables, the DSR problem 
can be re-formulated as: 

(^2) min - E Cn Z rrmSmn + r({£ mn }) 

(8a) 

subject to erg™ < er Gj „ < erg a *, Vra E S (8b) 

gn(<TG,n)+ E 4« <U = E ^ €n* (&) 



(L n Mt n Un < C, (m,») € £ (8d) 



feeAT„ 

rmax 



where Kirchhoff's current law ( |8c| ) is enforced at each node 
n E jV, and <j8bj are box constraints for the power supplied by 
controllable DG units. Matrices {Z mn } are positive definite, 
while matrices {M mn } are symmetric, positive semi definite, 
with rank 2. Thus, (P2) is convex, and it can be solved 
optimally via general-purpose interior point methods. What 
is more, (P2) can be reformulated as a second-order cone 
program (SOCP), and thus solved using efficient primal-dual 
interior point methods tailored to SOCP (see e.g., ED). 

The role of A in r({£ mn }) is to control the number of vec- 
tors {C ra J(m,n)e£ B (and, hence, of currents {w}(m,n)££ R ) 
that are set to zero. When A = 0, all branches £r are traversed 
by a non-zero current. Then, with A increasing, the number of 
lines where no current is flowing increases [14|. This implies 
that by adjusting A one can obtain either meshed topologies 
(low values of A), weakly-meshed, or even radial systems 
(high values of A). To rigorously show this, results from 
duality theory [22] are leveraged next to derive closed form 
expressions for the optimal line currents. 

Let {p n } and {p mn } denote the multipliers associated 
with ( [8c] i and ( |8d) , respectively, and consider the (partial) 
Lagrangian function of (P2), namely: 



1 x , 



T 

mn mnS mn 



, Z 7r ,n^ r . 



r ({Ll) 



A 



E^(gn(^,n)+ E4n%'»- E U 



n=l 

^ y y ' Pmn 

(n,m)e£ 4>ev mn 



Smn ?rmSr. 



(9) 



where £ := {£ mn }, &g '■— {vG.n}, and likewise p,p collect 
all the dual variables for brevity. Given ((9), the dual function 
and the dual problem take the form 



V{n,p) 



V opt = max V(a,p). 
{p n ho},p 



£($,a G ,ti,p) (10) 
(11) 



Since (P2) is convex, if there exists a feasible solution er G 
such that £LM m „£ m „ < C« for all (m,n) e £ and 



erg 1 " -< er G „ -< crg^ for all n E S (that is, Slater's condition 



£,°mn — argmin -j^l\ n Z mn ^ mrL + A||£ mn || 2 — Pmn^mn (12) 



holds), then (P2) has zero duality gap E2l Ch. 6]. Suppose 
that this is the case, and let C^^T and p opt ,p opt denote 
the optimal primal and dual solutions, respectively. The La- 
grangian optimality condition [22 Prop. 6.2.5] asserts that £ opt 
and er^* are also the minimizers of (|T0]> for p = p opt and p = 
p opt ; that is V op \p op \ p opt ) = £(£ op \ cr op \ p opt , p opt ). 
Thus, re-arranging terms of the Lagrangian function in a 
convenient way, and exploiting the decomposability of 
it turns out that the optimal currents flowing on the phases of 
line (to, n) are given as the solution of the sub-problem: 

1 

in 

L 

where A = for lines (to, n) E £\£r (whereas A > for all 
lines £r), and 

Z m n '■= Z m n + ^ t ptnn ^mn (13) 

4>ev mn 

ii — A( m )T . opt _ X(n)T . opt 

H'mn ■— -"-mn Pm -"-mn rn ■ \ L ^> 

Since Z mn is positive definite and ^^pt{ n pt M mn positive 
semidefinite, it follows that Z mn is positive definite and 
invertible. Thus, based on ( fl2] i, the optimal line currents are 
obtained next. 

Proposition 1: Per line (to, n) E £\£r, the optimal cur- 
rents are given by 



: opt 



mn Pmn 



(15) 



Proposition 2: If (to, n) E £r is a single-phase distribution 
line, then the optimal current ttT = {Wtn}, Wtn}V 
on phase cf> is given by the following soft-thresholding vector 
operation 

f <p.opt _ OlMmnllj ~ M+ 

^mn ~ / i i _ _ _i \ Pmn \^) 



($l{Zt n } 



<j>,opt 
Pmn \\Pmn\\2 



where [a]+ := max{0, a}. For lines (to, n) E £r that are 
two- or three-phase, the optimal vector of line currents is 
obtained via the following shrinkage and thresholding vector 
operation 

-l 



t opt _ 



r] opt I 



{HM m „ll 2 >A} [1 



i opt Z 



X 2 



(17) 



where rf pt E M + is the solution of the scalar optimization 
problem 

-l 

Pmn- (18) 



7)>0 Z 



A 2 . 



mn T ~^1-2\V m m\ 



Proof. Consider first the case of single-phase distribution 
lines, where ( |T2] > boils down to 

^ = argmin \ (5R{Z*„} + p*g*) Cntmn 

~ Pmn€mn + Ml^mnh ■ (19) 

The solver of ( fl9] l takes the form £ mn = z p mn f° r some 
scalar z > 0. In fact, among all possible £ mn with the 
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same ^2-norm, the Cauchy-Schwarz inequality implies that 
the maximizer of ^„4„„ is colinear with (and in the same 
direction of) fi mn . Thus, substituting £ TOn = zfJ, mn into ( fT9| ) 
yields the following problem in the scalar z: 

\ + ptT) Z 2 \\Hmnh 

-z||/x m „|| 2 + A|z| (20) 



z opt — are min 
6 z>a 2 



where X\z\ can be replaced by Xz since z > 0. The necessary 
and sufficient condition for z to minimize ( |23j > is [23, p. 92] 

||/i m J| 2 <A, if^* = 

»{^„ } + Ptln pt - 1 1 Mm„ 1 1 2 + A = 0, if z°p V 

(21) 

which is satisfied by 

z° pt = 3 i— : [||/*mnlla-A] + . (22) 

Relations (|T7^> — (fT~8^> can be proved along the lines of lfT31 . 
Specifically, (jT2j is first equivalently reformulated as the fol- 
lowing quadratic program with a second-order conic constraint 



van ^mn<±rnn 



mn>*7nn 



subject to 



-t. 



-< o 



(23) 



Next, derive the (concave) dual problem of ( |23) , which 
amounts to 



max -^(Mrn 
subject to ||xll2 — 1 



Ax)' ran (Prn It 



(24) 



where the constraint x G range(Z mn ) is left implicit, and x 
is the multiplier associated with the conic constraint in ( |23] l. 
Consider then the Lagrange dual of d24l, namely 



(25) 



min max -4/i mn + Ax) r Zj><>mn + Ax) 
r/>0 x ^ 



Recalling that Z mn is invertible, it turns out that the optimal 
solution of ( f2"5j l is given by 

-l 



2 ( 



A*. 



Mr, 



(26) 



with 77 as in (TTSj. Notice that the eigenvalues of (if pt Z mn + 
4^-1) denoted by {9i} for brevity, satisfy the inequality < 
6i < 2/A 2 for all i = 1,...,2\V„ 

A, the term (l/2)/x^„ ( 
is negative, thus implying that rf pt = 0. 



, n |. Thus, when ||/Lt, 



Mr, 



in ([18]) 

□ 



Some comments are now due in order to interpret the role 
of the multipliers and {ptf pt } in view of Ohm's Law, 

and to better appreciate the merits of the sparsity-promoting 
regularization term r({£ mn }). Notice first that from the com- 
plementary slackness condition El Prop. 6.2.5], one has that 
whenever the corresponding constraint d8dj> is not 



„<p,opt 



active. Suppose temporarily that pf,{° pt = for all lines, in 

I 2 m{Z mn } [cf. ©]. 



Since currents and voltages abide by Ohm's Law, (fl"5j)-(p"6l) 
imply that the legitimate unit for {fi mn } is the volt. In par- 
ticular, comparing ( p"5| ) with ([TJ reveals that jtx^ corresponds 
to the electrical potential difference between two nodes m 
and n connected by a line with a resistive matrix 5R{Z mTt }; 
that is, = [5R r {3fi{Z„ m }i m „},3 r {K{Z mn }i mn }] r . In 
other words, fi mn represents the contribution to the potential 
difference v m — v„ that is due to the resistive part of Z mn . 

With this connotation of fi mn , it follows from ( [T6| that 
single-phase lines equipped with a switch can be character- 
ized by a "virtual" line resistance, whose values is given 
by ^{ z tn) amplified by a factor ||jU ro J| 2 /[|| 
This resistive boost discourages high currents on line (m, n), 
something that in compressive sampling is usually referred 
to as "shrinkage operation" lfT4l . lfT5ll . Eventually, when 
II A*mn II 2 < X m e value of this virtual resistance goes to 
infinity, thus resembling an open switch. Notice further that 
the thresholding operator [||/- t mnll2 — A] , naturally suggests 
the order of magnitude of the parameter A that has to be 
used to (de)select a line. Although less intuitive, this shrinkage 
and thresholding operation effected through A||£ mn ||2 can be 
noticed also in ( [17] ) for lines with two and three phases. Here, 
the design variable 7/ is expressed in watt. Finally, variable 



<p,opt 



can be interpreted as an additional "virtual" resistance 
added to the conductor <f> of line (m, n) when the currents 
reaches its maximum allowable value /™^ x . In principle, since 
this extra resistive value introduce an additional virtual active 
power loss, higher values of lf nn are discouraged. 

The optimal topology of the network is obtained by dis- 
carding the distribution lines with an associated zero cur- 
rent. That is, £ opt := £\{(m,n) £ £ R : Cfn = °l- 
Further, the optimal line currents are readily computed as 



•.opt 



[£°mn]l:\V„ 



.1 



the voltage vector vi at the substation, and the optimal set 
of distribution lines £ opt , voltages {v„}^ =2 can be readily 
obtained by using ([TJ iteratively. Although (P2) is grounded 
on the affine approximation (|6j, it will be shown in Section [V] 
that the deviation between the load profile and the output 
of the DSR problem is on average less than 2 kW for the 
IEEE 37-node feeder fT6l - something on the order of the 
load forecasting error. Further, tests on this system revealed 
that it is possible to solve (P2) in less than 1 second (as 
opposed to several seconds or minutes required by solvers 
of mixed-integer nonlinear programs and genetic algorithms), 
thus rendering the proposed method particularly suitable for 
e.g., fast post-outage power delivery restoration purposes [ 17|. 

Finally, notice that in order to encourage the use of 
specific lines (m, n) 6 £r, the regularization function |7]) 
can be replaced by its weighted counterpart r'({£ mn }) '■= 



E 



(m,n)e£.R 



AmnHCr, 



with X mn > for all (m,n). For 



example, if the use of a line (j, k) is inadvisable due to e.g., 
ongoing maintenance or security concerns, a higher associated 
weight Xjk > X mn should be selected. 

IV. Accounting for voltage constraints 

Similar to various DSR renditions, the objective of (P2) is 
to obtain a topology that is likely to yield the lowest power 
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losses for a given predicted load profile. Based on the resultant 
optimal configuration {J\f,£ opt ), voltages and currents are 
fine-tuned in a subsequent optimization stage where more 
sophisticated techniques such as OPF are employed (see 
e.g., [24 1 and references therein). In some cases however, 
it may be desirable to introduce voltage regulation-related 
constraints in order to avoid network configurations that can 
potentially yield infeasible OPF solutions (meaning a set of 
voltages returned by the OPF solver not within prescribed 
minimum and maximum utilization limits). 

To effect voltage regulation, consider introducing a con- 
straint V$ £ B% per node n £ /V\{1} and phase <j) € V n , 
where Bf, is a given closed set collecting the admissible 
voltages (see e.g., [10], [19]); set B^ will be exemplified in 
the ensuing Sections |IV-A| and IV-B One way to enforce 
these constraints is to let the voltages become optimization 
variables, and formulate a joint DSR and OPF problem as in, 
e.g. iflOl . However, it is not convenient here to have voltages 
as optimization variables because: 

i) equality |5]) is nonconvex (as in OPF problems); and, 

ii) when binary variables are used to model the states of 
switches ||6l, flZl, ||8l, fl9], HOI, (TH, and constraints (|4cJ are 
employed 0, ED, flU, the solver would set i^* = for 
lines with x mn = 0. However, it is clear from ([TJ) that imposing 
' l mn = requires equating voltages at the two end points of 
line (m,n); that is, [v m ]-p mii ~ [ v n)v mn - This artifact renders 
the joint DSR and OPF problem infeasible in various practical 
cases. Consider for example a network with 5 nodes, lines 
£ = {(1,2), (2, 3), (2, 4), (2, 5), (3, 5), (4, 5)}, and switches 
in £r = {(2, 3), (2, 4), (2, 5)}. Suppose that only one switch 
must be closed in order to obtain a radial network, and this 
switch is the one on line (2,3). However, by setting v 2 = 
V4 = V5, the load demands at nodes 4 and 5 would not be 
satisfied, since no power is delivered at these nodes. 

One approach to resolving this issue is to discard (j4cj and 
replace {Z mn }( m ,„ )efH with {Z mn x mn }( min ) e£R . However, 
the resultant formulation yields a challenging bilinear problem 
with integer variables. Aiming at a computationally affordable 
DSR scheme, voltages are treated here as latent problem 
variables as shown next. 

Recall from Section |II-B that voltages at the substation 
V\ := pft^lvi}, Q^jvi}] 7 are typically taken as a reference 
for the phasorial representation lfl9l . Consider the network in 
Fig. |TJ Given v\, and assuming that the arc (1,2) originates 
at node 1 and ends at 2, voltages at node 2 can be expressed 
as [cf. {i}] 



3?{Z 12 } 
} 



3{Z 12 



} 

3?{Z 12 } 



-3{Z 12 



2 ■ 



(27) 



= * 12 



Likewise, if the sequence of nodes Wi^„ := {1, ... ,71} forms 
an undirected path 1 — >• 2 —> ... — >• n from the substation to 
node n, and none of the lines (m, n) : to, n € Wi^„ is 
equipped with switches, then v n can be expressed as v n = 



(m,n):m,nGWi_ 



W*mn£ m „, where a mn = 1 if the 



path traverses the directed edge (to, n) (which goes from m 
to n) in the opposite direction, and a mn = — 1 otherwise. A 



similar approach was taken in [6|. Based on these relations, 
voltage regulation can be readily effected by adding to (P2) 
the following constraint per node n: 



(m,n):m,n£Wi 



inn^jnn^mn ^ Bn 



(28) 



Once (P2) is solved, voltages are computed based on £ opt 



as shown in Section III When switches are densely deployed, 



there may not exist an undirected path connecting the substa- 
tion to a node n that includes only lines in £\£r. If this is 
the case, the substation must be replaced by another point of 
reference; that is, a node where the value (or an approximate 
value) of the voltages can be unequivocally determined. From 
the model set forth in Section |II-B| (see also lfT8l ). it can be 
noticed that approximate values of the voltages are readily 
available for nodes with a nonzero load demand; in fact, given 
the approximate injected currents u n in |6]), a reference voltage 
can be obtained via Q. Hereafter, the reference node will be 
generically denoted by n re f. Two possible choices for the set 
B n are presented in the ensuing subsections. 



A. Box constraints 

Let and VTjj be upper and lower bounds, respectively, 
on the real part of voltage Kf. Likewise, let V„c* and VJc, 
denote the counterparts for GjV^f}. Then, upon collecting 
these quantities for all phases <f> £ V n in the vectors i> n := 

[Ry^yr and *>„ : = i^y ^yr. m «° 

be readily reformulated as follows: 

(P2') min - V iXin^mnimn + r ({£rwJ) 



subject to ( |8b| t, ([8c]), (|8d]i, and 

(m,n):m,nGW T i T , e y 



Problem (P2 r ) is convex, and thus efficiently (and optimally) 
solved via either general-purpose interior point methods, or, 
primal-dual schemes tailored to SOCP ||2TI . Through ( |29] > it 
is possible to constrain both the magnitude of V* and its 
deviation from the nominal phase ZY^ . Further, introducing 
these additional constraints does not alter the expressions for 
the optimal line currents provided in Propositions [TJ and [2] 



B. Nonconvex constraints on voltage magnitudes 

In traditional OPF approaches, it is usual to consider lower 
and upper bounds on the voltage magnitudes, which habitually 
coincide with minimum and maximum utilization and service 
voltage levels [19|. Consider then introducing in (P2) the 
constraints V n < \V$\ < V n , per node n and phase <fr £ V n . 
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This leads to the following DSR reformulation: 

(-P2") m )n 2 E + r({€ m „}) 

subject to < |8bl >, ( f8c] ) ? ( |8d] >, and 

r 



\ m.n / 



; m n ^ mnsmn 



<V*,<j)£ V„ (30) 



where := I2 (8 diag(e^). Unfortunately, the constraint 

- (i/„ re/ + ^£,) T Qt (v„ nf + ¥£) +f„ 2 < 0, where * and 
£ are defined in the obvious way, is nonconvex since the 
function on the left hand side is concave. However, this source 
of non-convexity can be efficiently addressed by resorting to 
convex approximation techniques. Among candidate methods, 
the successive convex approximation (SCA) approach pro- 
posed in [25] is well suited for the problem at hand because it 
guarantees first-order Karush-Kuhn-Tucker (KKT) optimality 
under mild regularity conditions. 

To briefly illustrate the general SCA method, consider an 
optimization problem 

min/o(£) (31a) 
subject to < 0, k = 1, 2, . . . ,K (31b) 

where /o(£) is convex and differentiable, k = 1, . . . , K, 

are differentiable functions, and the feasible region T := {£ G 
^|/fc(£) < 0, k = 1, . . . , if} is compact. Then, starting from 
a feasible point £' ' 6 T , a series of approximate problems 
can be solved to locate a KKT optimal point of the original 
(non-convex) problem. For each k = 1, . . . , K , let C ) 
denote the surrogate function for which may depend 

on the solution to the problem in the (j — l)-st iteration. 
The approximate problem to solve in iteration j is 

min/o(0 (32a) 

subjectto/ fc (C;C W) )<0, k = l,2,...,K (32b) 

whose feasible region is denoted as J 7 ^'. Provided that 
) satisfies the following conditions cl)-c3) for each 
k = 1, . . . , K, the series of solutions j — 1,2,.. ., to 
the approximate problems converges to the KKT point of the 
original problem ( (31) : 

ci) /*«)</*(€; V^J« 

c2) = /fc^C^) 

c3) V/fe(^)=V/fe(^;^). 

In order to apply the SCA method to (P2"), an 
appropriate surrogate constraint for the nonconvex 
lower bound in |30) needs to be determined. Since, 

- {y n ^s + *l) T Qt ("n ref + *£) is quadratic and 
concave, a linear approximation around the feasible point 
{£m«} satisfying cl)-c3) can be readily found. In fact, after 
standard manipulations, it follows that the convex problem to 



be solved at iteration j of the SCA algorithm is: 

(P2"W) min \ £ CnZmntmn + ratmn}) 

a,n} I 

(n,m)Et 

subject to (|8b]), ([8c|,([8d]), and 

(*£(Z U) )) T {e U) - - + V% < (33a) 



*^) r Q*K. o/ + *0 <Vl (33b) 



with 



a£(£«) := 2* T Qt (V„,. c/ + (33c) 

:= + ^ (J) ) T Q^ (iV., + • (33d) 

The resultant SCA algorithm amounts to solving a sequence of 
convex problems (P2"( )) ... -> (P2"W)) -> where the 
feasible point is taken to be the solution of (P2"^ _1 ^) 
(and £(°> is an arbitrary initial feasible point). Convergence of 
this scheme is formalized next (see also GBI ). 

Proposition 3: Given an initial feasible point the iter- 
ates generated by (P2' ,(j )), j = 0, 1,2, . . ., converge to a 
KKT solution of the nonconvex DSR problem (P2"). 

In the ensuing section, the proposed DSR algorithms are 
tested on the IEEE 37-node feeder fl6l . 

V. Numerical Experiments 
Consider the IEEE 37-node feeder [16], which is a portion 
of the three-phase 4.8 kV distribution network located in 
California. Compared to the scheme in [16|, 8 additional 
three-phase lines equipped with sectionalizing switches are 
considered, as shown in Fig [Tj The parameters of these 
additional lines are listed in Table [I] where the admittance 
matrices corresponding to the configuration indexes 723 and 
724 can be found in |16|. The line impedance matrices for the 
original lines are computed as specified in lfl6l . A balanced 
load of 85 kW and 40 kVAr per phase is added to node 
23 to represent an additional residential demand. Further, 
DG units are located at nodes 9, 13, 15, 19, 26, 32, and 36. 
These units operate at unity power factor, and can supply a 
maximum power of 50 kW per phase. Two different setups 
are considered, depending on the number and locations of 
switches: 

1) Test 1: E R ,i = {(8, 14), (6, 20), (10, 16), (20, 26), 
(16, 24), (10, 17), (24, 33), (26, 35)}; and, 

2) Test 2: £ R>2 = £ R>1 U {(17, 22), (23, 24), (23, 25), 
(29,30)}. 

The voltage at the distribution substation is set to Vi = 
[1Z0°,1Z - 120°,lZ120°] r pu. Finally, box constraints on 
the voltages are considered; specifically, the lower and upper 
bounds of the real and imaginary parts of the voltages are 
such that K{V^} G - 0.0354, lSt{Vf] + 0.0354] pu 

and %{V£} G [%{V{} - 0.0354, 3{F/} + 0.0354] pu. This 
translates to having the magnitude of the voltages in the range 
I V* I G [0.95,1.05] pu. 

The optimization package CV>|^] along with the interior- 
point based solver SeDuMij^Jare employed to solve the DSR 

2 [Online] Available: http://cvxr.com/cvx/ 

3 [Online] Available: http://sedumi.ie.lehigh.edu/ 



8 




Fig. 1. Modified IEEE 37-bus test feeder. 



problems in MAT LAB. The average computational time re- 
quired by SeDuMi to solve the DSR problem was 0.4 seconds 
on a machine with Intel Core i7-2600 CPU @ 3.40GHz. 

TABLE I 

Additional lines in the modified IEEE 37-node feeder 



Line 


Conf. 


Length (ft) 


Line 


Conf. 


Length (ft) 


(8,14) 


723 


1144 


(16,24) 


724 


1580 


(6,20) 


724 


1320 


(10,17) 


724 


1137 


(10,16) 


724 


847 


(24,33) 


724 


1315 


(20,26) 


724 


815 


(26,35) 


724 


377 



Fig. [2] depicts the magnitude of currents flowing on lines 
equipped with switches, for different values of the sparsity- 
promoting parameter A. Three rows per line are reported, 
where the first row corresponds to phase "a," and the third one 
to phase "c" (all lines are three-phase). The voltage magnitude 
is color-coded, where white represents a zero current (that is, 
an open switch), while red hues are used to capture different 
values for \I%\ > 0, in pu. In the upper plot (a), it can be 
clearly seen that the number of open switches increases as 
A increases, and the original tree topology described in |[T6ll 
is obtained at saturation. Interestingly, the first lines that are 
discarded are (10, 16) and (10, 17), which implies that the 
majority of the power supplied to that part of the network 
comes from the DG units. As expected, similar trends are 
evidenced in Fig. [2jb) for Test 2, thus further corroborating 
the merits of the proposed method based on sparsity-promoting 
regularizations. It can be noticed that at saturation only four 
switches are left closed and, interestingly, the resultant tree 
topology is different than the original one in fl6l . In fact, the 
switches in (23,24), (23,25) are open, while lines (16,24), 
(20, 26), and (24, 33) are used to deliver power to the loads. 

Fig. |3]reports the expected active power loss as a function of 
A, for both cases. Notice first that power losses are in general 
lower in Test 2. Then, it can be clearly seen that the power 
loss monotonically increases as A increases (that is, with the 
number of open switches). This motivates augmenting the cost 
of (P2) with a term that accounts for the possible maintenance 



(6,20) 
(8.14) 
(16,24) 
(10,16) 
(10.17) 
(20,26) 
(24,33) 
(26,35) 



icj [p<«] 



20 25 30 35 



(a) Test 1. 



45 50 



lu 0:5 

1 0.045 
I 0.04 
I 0.035 

-] 0.03 
1 0.025 



(6,20) 
(8,14) 
(16,24) 
(10,16) 
(10.17) 
(20,26) 
(24,33) 
(26,35) 
(17,22) 
(23,24) 
(23,25) 
(29,30) 



(b) Test 2. 

Fig. 2. Absolute value of the currents on lines Sr, for different A. 

costs of lines and switches, in order to find a possible trade- 
off between active power loss and number of utilized lines. If 
this economical cost incurred by the use of a line is modeled 
by using a strictly convex function, the optimal A can be 
efficiently found via bisection. 

Since the linear relation ^ introduces an approxima- 
tion error when PQ loads are present, the average devi- 
ation from the nominal loads is quantified next. Specifi- 
cally, the deviation for the real power is defined as AP = 
(V E„ \K I) En E \P£~H\> where Vn}s a set collecting 
the phases at node n with a non-zero load; P£ is the output 



of the DSR scheme; and P$ = P 



G.n 



L,1 



is the true 



real power that would be obtained by considering the exact 
nonlinear relation Q. The counterpart AQ is defined in a 
similar way. These deviations yield only AP = 1.49 kW 
and AQ = 0.92 kVAr for Test 1, and AP = 1.35 kW and 
AQ = 0.85 kVAr for Test 2. The error is very small, and does 
not jeopardize optimality of the obtained topology (as would 
happen in case of higher approximation errors). 

The magnitude of the voltages V r f are reported in Fig. [4] 
for different values of A. The substation voltages are not 
included. Voltages are computed as explained in Section III, 
once the optimal topology has been identified. As expected, 
the "cloud" formed by all the voltage mag nitudes {\V*\} 
is confined within the considered lower and upper bounds, 
which are marked with red dotted lines. Fig. [5] reports the 
active power supplied by the DG units. It can be seen that in 
both cases the majority of DG units supplies 50 kW, which 
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Fig. 3. Overall active power loss [kW]. 
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(b) 

Fig. 4. Voltage magnitude at the nodes [pu]. (a) Test 1. (b) Test 2. 

corresponds to the maximum power that can be generated. In 
Fig. [5ja) the power supplied by DG units 7 and 8 is smaller, or 
even zero, on phase "b." This is because the load of that part 
of the network on phase "b" is relatively low, and using those 
DG units to supply loads that are electrically far from nodes 
32 and 36 incurs higher power losses on the lines. Similar 
conclusions can be drawn for Fig. |5jb). However, in this case, 
because of the new topology, DG units 7 and 8 are judiciously 
utilized to supply other loads. 

VI. Concluding Remarks 

A DSR problem was considered for three-phase distribution 
systems featuring DG. Leveraging the notion of group-sparsity, 
and adopting an approximate linear relation between powers 
and injected currents, a novel convex DSR formulation was 
proposed. Being convex, the proposed DSR problem can be 
solved efficiently even for distribution networks of large size. 
The ability of the proposed scheme to select the topologies 



Fig. 5. Active power supplied by the DG units [kW]. (a) Test 1. (b) Test 2. 

that minimize the overall active power loss was demonstrated 
via numerical tests, and it was also justified analytically. 

Appendix 

Fig. [6] reports magnitude (in pu) and phase of the injected 
(absorbed) currents for the IEEE 37-node feeder. The orig- 
inal scheme described in |[T6l is considered, with voltages 
and currents obtained by using the OPF solver proposed 
in |24l . Black, green, and orange circles correspond to the 
exact values of the currents for phases a, b, and c, respec- 
tively. These values of current magnitudes and phases are 
compared with the ones obtained by using the approximate 
relation (|6l, which are marked with "x". It can be seen that 
the incurred approximation error is very small (see also the 
results reported in lfT8l ). More precisely, the average error 
A := (l/EJ^DEnE^f - It\, where /£ m is the 
current computed using |6]l amounts to 0.1767 Ampere. 

This approximation error is certainly acceptable for two 
reasons: first, using |6} in conjunction with compressive sam- 
pling methods, gives rise to the formulation of a convex DSR 
problem for a system featuring PQ loads. Secondly, voltages 
and currents are usually adjusted in a subsequent optimization 
stage via OPF. 
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